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Abstract 

An idealized multigrid algorithm for the computation of propagators of staggered 
fermions is investigated. Exemplified in four-dimensional SU{2) gauge fields, it is shown 
that the idealized algorithm preserves criticality under coarsening. The same is not true 
when the coarse grid operator is defined by the Galerkin prescription. Relaxation times 
in computations of propagators are small, and critical slowing is strongly reduced (or 
eliminated) in the idealized algorithm. Unfortunately, this algorithm is not practical for 
production runs, but the investigations presented here answer important questions of 
principle. 
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1 Introduction 



In Monte Carlo simulations of lattice gauge theories with fermions the most time-consuming 
part is the computation of the gauge field dependent fermion propagators. Great hopes to 
compute propagators without any critical slowing down (CSD) are attached to multigrid (MG) 
methods [||, 0, |, ^, 0, ||, |ll|. However, up to now no practical MG algorithm has 

been found for fermions. 

In this note an idealized MG algorithm is investigated for staggered fermions in four- 
dimensional SU (2) gauge fields. It will be shown that the idealized algorithm preserves critical- 
ity under coarsening, which is not true when the coarse grid operator is defined by the Galerkin 
prescription. This finding explains the failure of simple variational-like MG methods, at least 
for algorithms with nonoverlapping blocks or trivially overlapping blocks. 

Relaxation times in computations of propagators with the idealized MG algorithm are small, 
and CSD is strongly reduced or eliminated. Unfortunately, this algorithm is not practical for 
production runs, but the investigations presented here answer important questions of principle. 



2 Multigrid Method 

For given / we wish to solve an equation 

DoX = f with Do={-P' + m') (1) 

by MG, where ^ is the gauge covariant staggered Dirac operator, and m is a small quark mass. 

The following MG notations will be used. The fundamental lattice is denoted by A°. The 
first block lattice is obtained by coarsening with a factor of L^. Thus has fewer 
sites than A° (in d space-time dimensions). Restriction and interpolation operators C and A, 
respectively, are given by kernels C{x,z) and A{z,x) with z G A'', x G A^. Note that C{x,z) 
and A {z, x) are Nc x Nc matrices in a gauge theory with N^. colors. Also, C and A depend on 
the gauge field, although this is not indicated explicitly. 

We use a blocking procedure for staggered fermions which is consistent with the lattice 
symmetries of free fermions 0. This forces us to choose Lb = 3. Even Lb are not allowed. In 
four dimensions, coarsening by a factor of three reduces the number of points by 81. Therefore 
only a two-grid algorithm was implemented. The residual equation on the coarse grid was 
solved exactly by the conjugate gradient algorithm. 



The averaging kernel C is chosen according to the ground-state projection definition [|1^, 
13, H, H. In the present work C fulfills the gauge covariant eigenvalue equation(s) 



{-A^,,C*){z,x) = \oix)C*{z,x) (2) 

together with a normalization condition CC* = 11, and a covariance condition C{x,x) oc 11 
where x denotes the center of block x. In Eq. (|]), Xo{x) is the lowest eigenvalue of — Atv.x, 
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and — Atv.x is the gauge covariant fermionic "two-link lattice Laplacian" - defined through 
ip^ = A + a^iyF^„ - with "Neumann boundary conditions (b.c.)". Neumann b.c. means that 
terms in A are omitted where one site is in block x and the other one is in a neighboring block. 

The ground-state projection method is numerically implementable in four-dimensional non- 
abelian gauge fields [Q, and since the method is gauge covariant, no gauge fixing in com- 
putations of propagators is required. For staggered fermions in non-abelian gauge fields two 
qualitatively different proposals were made for ground-state projection 0. We call these pro- 
posals "the Laplacian choice" and "the Diracian choice". The Laplacian choice is the one 
described above. In the Diracian choice one substitutes a block-local approximation of for 
A in (0). This latter choice would be superior because it takes also the field strength term 
Ffj,iy into consideration. However, it was proved numerically P, |l^ that the Laplacian choice 
for C defines a good blockspin in arbitrarily disordered gauge fields.^ For this reason only the 
Laplacian choice of C has been implemented yet. 



3 Spectrum of -^^ and CSD 

The square of the staggered Dirac operator (plus mass term) couples only even lattice sites to 
even sites, and odd sites to odd sites. Therefore the matrix elements of — can be arranged 
in such a way that —]/)^ can be written symbolically as 



(3) 



where —-^even/odd — restricted to the even/odd sublattice. 

Let us denote the spectrum of — by S (— ^^). It equals the union of the spectra of — -^even 
and S {—I/)'^) = S {—^l^^j^)\JS {—^l^^). The spectra are gauge invariant. Moreover, 

for any gauge field configuration one has the equality 

S i-PlJ = S (-^odd) ■ (4) 

A simple proof of is as follows. Consider the lattice operator —-^even/odd a block matrix 
with Nc X Nc elements — -^even/odd (-^i) -^2)- The matrices — -^even ~-^odd similar, and 
therefore they have the same spectrum. We recall that two matrices A and B are called similar 
if there exists an invertible matrix T such that B = TAT~^. (We also recall that if v is an 
eigenvector of A with eigenvalue A, then Tv is an eigenvector of B with eigenvalue A.) In the 
case considered here, we look for a lattice operator T with the property 

5:T(t.,z')(-^e\en)(^',^) = E(-^odd) («^, «^')T(«;', ^) (5) 



^■'For every averaging kernel C there exists an associated ideal interpolation kernel A] see Sec. ^. C defines 
a good blockspin if this A decays exponentially. 
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where 2;, z' and w, w' denote even and odd lattice sites, respectively. Eq. (^) is fulfilled if we 
choose the matrix elements of T to be T(iy, z) = Jp{w^z). This choice of T is not invertible in 
pure gauges, but in that case the equality (H) of spectra is obvious anyhow. 



In Ref. U5[ a more complicated proof of (^) was given which uses an analyticity argument 
in connection with a hopping expansion of (— + m^)~^ for large mass. 

In Refs. |TU|, ^ the author pointed out that in conventional relaxation algorithms for 
propagators of staggered fermions there exists a scaling law for relaxation times r which reads 

T = — — - With Am = m — m^^ (6) 

for small Am^, where m^^. is the lowest eigenvalue of and const, is independent of the 

lattice size. For bosonic propagators the validity of is known analytically and has also 
been confirmed to a high accuracy numerically |0, |ll| . 

A consequence of (§) is that in conventional relaxation algorithms for staggered fermions, 
CSD will be the same on the even and the odd sublattice. 



4 Idealized Multigrid Algorithm 

Up to now no practical MG method has been found for fermions. Mack pointed out that it 
is essential for fighting CSD in MG computations that interpolation kernels should be smooth 



Iq . This requirement is not fulfilled in MG algorithms where one uses gauge covariant gener- 
alizations of piecewise constant interpolation with nonoverlapping blocks. Mack suggested an 
interpolation kernel ,4 as a starting point for numerical work which was used successfully by 
Gaw§dzki and Kupiainen in constructive quantum field theory An idealized MG algorithm 
using the natural gauge covariant generalization of the Gaw§dzki-Kupiainen kernel had been 



investigated numerically in four-dimensional SU (2) gauge fields for bosonic propagators []^, [TT 
There CSD could be eliminated completely. Here we turn to an idealized MG algorithm for 
staggered fermions. This algorithm will not be practical for production runs, but it is important 
to answer questions of principle, and to recognize the features which a successful method must 
have. 

Given the averaging kernel C, there exists an ideal choice of the interpolation kernel A. It is 
determined as follows. For every function ( "block spin" ) $ on A^, = ^ $ minimizes the action 
< (j) , Do(j) > subject to the constraint C(j) = For the purpose of numerical computations, it 
is convenient to determine the optimal A as the solution of the equation 

{[-^^ + + kC*C]A) {z,x) =kC*{z,x) (7) 

for large k. C* denotes the adjoint of C. The layers of an MG decouple completely when this 
A is used for interpolation, and when coarse grid operators Di are defined as C{—jp'^ +m^)^. 
These coarse grid operators are automatically hermitean and equal A*{^—JP^ + rr?^A . 
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11 . With the definition of 



The optimal A is favored by an argument of dynamics 
smoothness that covariant derivatives are small, the above characterization of A as solution 
of an extremization problem can be rephrased: that A is the smoothest interpolation kernel, 
subject to the constraint CA = 1. 



5 Numerical Results 

Because of the considerations of Sec. || we decided to investigate the ideahzed two-grid algorithm 
only on the even sublattice. This restriction mitigates also the storage space requirements for 
the ideal A a little bit. Note that Eq. can also be broken up into an equation for the 
interpolation kernel on the even sublattice and one on the odd one, if the averaging kernel C 
does not mix even and odd sites. This requirement is fulfilled both for the Laplacian and for 
the Diracian choice of Ref. 0. In both proposals the coarse grid sites can be separated into 
even and odd sites, and the ideal effective Diracian —C^^A (as well as the Galerkin operator 
—C^'^C*) can be decomposed analogously to (R|). 



Here we made the Laplacian choice and computed C by the efficient algorithm of Ref. pA 
Numerical work was done in SU{2) lattice gauge fields on 6^ and 12^ lattices, covering all 
possible values of /3 = A/g^ between oo and zero. The system (|^) was solved by means of the 
conjugate gradient algorithm where iterating was stopped when the RMS residual was less than 
]^g-io_ rpj^g statement that —CJp^A is automatically hermitean was confirmed up to round-off 
errors of order 10~^ or less. 



5.1 Lowest Eigenvalues 

Let us first look at the lowest eigenvalues of —1/)^ on the fundamental lattice A'', and see how 
they are transferred to the block lattice A^. The role of this transfer for the performance of the 
parallel-transported MG method of Ben-Av, Brandt and Solomon was pointed out in Ref. . 

First the lowest eigenvalues —ml^ of —]/)^ were determined by inverse iteration. This method 
allows to determine — wi^j. to an accuracy of 10"'' or better 0. Then optimal interpolation 
kernels A were computed as solutions of Eq. (0) with rn? = m^j. , and for k = 10^. Results 
for — m^j, and for the lowest eigenvalues of the ideal coarse grid operator C{—jp^ + rn^^-^A are 
given in Table |l[ The last column of Table |l| contains results for the Galerkin definition of the 
coarse grid operator. This operator is used in variational MG where interpolation is done by 
C* . (The Galerkin operator retains the locality properties of —jp^ in arbitrary gauge fields.) 

One sees that for any value of the gauge couphng the idealized algorithm maps a critical 
system on A'^ onto a critical system on the block lattice. In contrast, variational MG does 
not have this property. The Galerkin operator C{—Jp'^ + m^j.)C* is far from being critical in 
nontrivial gauge fields. 

The results of Table |1| supply another explanation for the failure of variational MG which 
was ascertained in Refs. f^, [1^, |T^. One cannot expect that a (nearly) critical problem on A° 
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can be solved by means of an auxiliary problem with fewer degrees of freedom on A^, if the 
auxiliary problem is not critical as well. The Galerkin operator is only critical in trivial gauge 
fields,^ and only there is CSD eliminated by the variational MG method in computations of 
propagators. 

The effect of adding a mass term Am^ to (— + m^j.) is as follows. The values given in the 
last two columns of Table ^ are shifted by the amount determined by Am^. This is obvious in 
case of the Galerkin operator because the averaging kernel is normalized such that CC* = 1. 
In case of the idealized algorithm, CA tends to 1 for k — oo. For finite n one finds deviations 
from 11. Examples are given in Table |^. These deviations are small enough to have no effect 
in practice. When one computes the lowest eigenvalues of —C^^A, one recovers the negative 
critical masses of the gauge fields to the same accuracy as they are given in Table 0. Hence, 
the small negative values partly found for the lowest eigenvalues of C(— + mlj.)A are really 
due to numerical inaccuracies. 



5.2 Performance of the idealized algorithm 

Finally, we report results of computations of propagators by means of an idealized algorithm. 
We used coarse grid operators C{—jp'^ + m^)A with masses = m^j. + Am^, Am^ > and 
small. For all values of Am^ only one interpolation kernel A was used, namely the "critical" 
one which solves (0) with = m^^,. Actually, one should use an m^-dependent ^-kernel, viz. 
the solution of with being the mass under consideration. However, in case of bosonic 
propagators the procedure described here was successfull 0, and therefore we used it also 
first attempt in case of staggered fermions. 

Tables ^ and ^ comprise results for relaxation times as a function of Am^ on 6^ and 12^ 
lattices, respectively. The values given for the relaxation parameter u are optimal within ±0.05. 
(For the configuration on the 12^ lattice at (3 = 2.5, uj = 1.72 is close to optimum, while the 
results for u = 1.65 are given as additional information.) 

We rediscover here an observation which was made earlier [0, 11 1: The use of a relaxation 
parameter u different from one in MG computations contradicts the conventional wisdom. 
According to this wisdom the only job of the relaxation procedure on A° is to smoothen the error, 
and this job is well done by Gauss-Seidel iteration. The conventional wisdom was confirmed 
by numerical results in trivial gauge fields [0. However, the picture changes for propagators 
in nontrivial gauge fields. 

In order to be sure about the correct determination of the value of m^^,, it was checked that 
conventional SOR and the variational MG algorithm (with the Laplacian choice of C) both 
exhibited CSD, i.e. relaxation times r follow perfectly the scaling law (^. const, in (^ is of 
order one. 

The results of Tables ^ and § show that the 1/Am^ divergence of r's on lattices of a fixed 

^^In exact arithmetics the entries for f3 = oo (reahzed as random pure gauge fields) in Table |l| would be zero 
for any lattice size. 
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size is eliminated in the idealized MG algorithm. Relaxation times are bounded and small. 
(Only for very disordered gauge fields at physically uninteresting values of the gauge coupling, 
r's are not so small, but nevertheless bounded.) We conclude from Tables |^ and ^ that CSD in 
computation of propagators is strongly reduced. It is hard to judge a possibly remaining volume 
effect, but one might be tempted to say that CSD can be eliminated for practical purposes, in 
principle. At this point one should also recall that in the investigations reported here, only one 
^-kernel was used for all values of m^. If one used an m^-dependent kernel, the results would 
probably improve further, at least they cannot become worse. 

Acknowledgments 

I wish to thank Gerhard Mack for many stimulating discussions. Financial support by 
Deutsche Forschungsgemeinschaft is gratefully acknowledged. The computations reported here 
were performed on the CRAY Y-MP of HLRZ Jiilich. 

References 

[1] R. Ben-Av, A. Brandt, and S. Solomon, Nucl. Phys. B329, 193 (1990); R. Ben-Av, A. 
Brandt, M. Harmatz, E. Katznelson, P.G. Lauwers, S. Solomon, and K. Wolowesky, Phys. 
Lett. B253, 185 (1991); Nucl. Phys. B (Proc. Suppl.) 20, 102 (1991); R. Ben-Av, P.G. 
Lauwers, and S. Solomon, Nucl. Phys. B374, 249 (1992); P.G. Lauwers and S. Solo- 
mon, Int. J. Mod. Phys. C3, 149 (1992); P.G. Lauwers and T. Wittlich, Inversion of 
the Fermion Matrix in Lattice QCD by Means of Parallel- Transported Multigrid (PTMG), 
preprint BONN 92-29 (September 1992), to appear in Int. J. Mod. Phys. C4; Nucl. Phys. 
B (Proc. Suppl.) 30 (1993) 261. 

[2] M. Harmatz, P. Lauwers, S. Solomon, and T. Wittlich, Nucl. Phys. B (Proc. Suppl.) 30 
(1993) 192. 

[3] R.C. Brower, C. Rebbi, and E. Vicari, Phys. Rev. D43, 1965 (1991). 

[4] R.C. Brower, C. Rebbi, and E. Vicari, Phys. Rev. Lett. 66, 1263 (1991); R.C. Brower, 
K.J.M. Moriarty, C. Rebbi, and E. Vicari, Nucl. Phys. B (Proc. Suppl.) 20, 89 (1991); 
Phys. Rev. D43, 1974 (1991); R.C. Brower, R.G. Edwards, C. Rebbi, and E. Vicari, Nucl. 
Phys. B366, 689 (1991). 

[5] A. Hulsebos, J. Smit, and J.C. Vink, Nucl. Phys. B (Proc. Suppl.) 20, 94 (1991); Int. J. 
Mod. Phys. C3, 161 (1992); Nucl. Phys. B368, 379 (1992); J.C. Vink, Phys. Lett. B272, 
81 (1991); Nucl. Phys. B (Proc. Suppl.) 26, 607 (1992). 

[6] T. Kalkreuter, G. Mack, and M. Speh, Int. J. Mod. Phys. C3, 121 (1992). 

[7] T. Kalkreuter, Phys. Lett. B276, 485 (1992). 



7 



[8] G. Mack, T. Kalkreuter, G. Palma, and M. Speh, in Computational Methods in Field 
Theory, Lecture Notes in Physics 409, ed. H. Gausterer and C.B. Lang (Springer Verlag, 
Berlin, 1992); M. Baker, T. Kalkreuter, G. Mack, and M. Speh, Neural Multigrid for 
Gauge Theories and Other Disordered Systems, preprint DESY 92-126 (September 1992), 
to appear in the proceedings of the 4th International Conference on Computational Physics 
PC92, Prague, August 1992. 

[9] T. Kalkreuter, Int. J. Mod. Phys. C3, 1323 (1993). 

[10] T. Kalkreuter, Nucl. Phys. B (Proc. Suppl.) 30 (1993) 257. 

[11] T. Kalkreuter, Multigrid Methods for the Computation of Propagators in Cauge Fields, 
Ph. D. thesis and preprint DESY 92-158 (November 1992), shortened version submitted 
to Int. J. Mod. Phys. C. 

[12] G. Mack, unpublished. 

[13] A. Hulsebos, J. Smit, and J.C. Vink, Nucl. Phys. B (Proc. Suppl.) 9, 512 (1989). 

[14] T. Kalkreuter, Nucl. Phys. B376, 637 (1992). 

[15] O. Martin, and S.W. Otto, Phys. Rev. D31, 435 (1985). 

[16] G. Mack, in Nonperturbative Quantum Field Theory, ed. G. 't Hooft et al. (Plenum Press, 
New York, 1988). 

[17] K. Gaw§dzki and A. Kupiainen, Commun. Math. Phys. 77, 31 (1980); ib. 99, 197 (1985). 



8 



Tables 



Table 1: Lowest eigenvalues of the negative squared Dirac operator for staggered fermions on a 
fundamental lattice A°, and the lowest eigenvalues of the ideal block operator C(— + mlj.)A 



of the Galerkin operator C(— + mlj.)C*. 






— m^j. = lowest eigenvalue of 


lowest eigenvalue of 


lowest eigenvalue of 




|A°I 




C{-P' + mi) A 


C{-P' + mi)C* 


oo 


6^ 


7.98-10-29 


3.09-10-21 


9.04-10-1^ 


5.0 


6^ 


0.0013413 


-1.40-10-11 


1.0341375 


3.0 


6^ 


0.3497739 


6.68-10-12 


2.8529730 


2.8 


6^ 


0.2441995 


7.27-10-11 


2.8412484 


2.8 


6^ 


0.2748178 


4.33-10-12 


3.2338422 


2.7 


6^ 


0.2004647 


-6.42-10-11 


3.4374150 


2.6 


6^ 


0.1740946 


-2.78-10-11 


3.1884409 


2.5 


6^ 


0.0698942 


-5.96-10-1^ 


3.0198011 


2.4 


6^ 


0.0010729 


-1.82-10-11 


3.1265394 


2.2 


6^ 


0.0007099 


9.50-10-11 


3.4488349 


2.0 


6^ 


0.0000732 


6.08-10-1° 


3.5620285 


0.0 


6^ 


0.0000287 


-1.21-10-1° 


3.7354123 


oo 


12^ 


8.15-10-29 


-2.33-10-1^ 


9.13-10-1^ 


3.0 


12^ 


0.0779810 


-7.96-10-12 


2.3190549 


2.7 


12^ 


0.0368447 


-1.15-10-11 


2.3967694 


2.5 


12^ 


0.0005742 


-1.03-10-1° 


2.6881098 


2.4 


12^ 


0.0001865 


-1.10-10-^ 


2.8248637 
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Table 2: Accuracy of CA = 11 for staggered fermions. \\CA — ]l||oo denotes the maximal trace 
norm of CA{x,y) — 6{x — y) over all pairs {x,y) of block lattice sites, and \\CA — 1\\2 is the 

RMS of these norms. 

(3 |A0| \\CA ~lh \\CA -l\\oo 



ooFI 

I—I 


6^ 


3.98- 


10- 


-15 


2, 


.13' 


■10" 


-14 


5.0 


6^ 


6.62- 


10- 


-7 


1, 


.87. 


•10- 


-6 


3.0 


6^ 


2.77- 


10- 


-6 


9, 


.43' 


•10- 


-6 


2.8 


6^ 


2.59- 


10- 


-6 


8, 


.66- 


•10- 


-6 


9 8 





O.Oy' 


iU 


-6 


1 

i , 


.U ( • 


• iU 


-5 


2.7 


6^ 


3.76- 


10- 


-6 


1, 


.18' 


•10- 


-5 


2.6 


6^ 


3.18- 


10- 


-6 


L 


.02. 


•10- 


-5 


2.5 


6^ 


2.43- 


10- 


-6 


8, 


.15' 


•10- 


-6 


2.4 


6^ 


2.39- 


10- 


-6 


6, 


.76- 


•10- 


-6 


2.2 


6^ 


1.51- 


10- 


-6 


4, 


A8- 


•10- 


-6 


2.0 


6^ 


9.40- 


10- 


-7 


2, 


.95- 


•10- 


-6 


0.0 


6^ 


2.84- 


10- 


-7 


9, 


.23- 


•10- 


-7 


oo 


12^ 


2.32- 


10- 


-6 


2, 


.32- 


•10- 


-5 


3.0 


12^ 


1.02- 


10- 


-6 


1, 


.\2- 


•10- 


-5 


2.7 


12^ 


1.00- 


10- 


-6 


1, 


W- 


•10- 


-5 


2.5 


12^ 


9.08- 


10- 


-7 


L 


.05' 


•10- 


-5 


2.4 


12^ 


7.92- 


10- 


-7 


9, 


.12. 


•10- 


-5 



"In case of the pure gauge on the 6^ lattice (and only in this case!), the ideal MG scheme is identical to 
the Galerkin definition with A — C* ("covariant piecewise constant" interpolation); since C is normalized as 
CC* = 1, the finite norms of CA — II in this case are due to round-off errors. 



10 



Table 3: Dependence of relaxation times r on Am^ in the ideal MG algorithm with lexicographic 
SOR for computing propagators of staggered fermions on 6^ lattices. (The two configurations 
at /9 = 2.8 are different.) 

r for Am^ = 



/7 




1 










10""' 


10"*^ 


oo 


1.17 


0.8 


0.8 


0.7 


0.7 


0.7 


0.7 


0.7 


5.0 


1.35 


1.2 


1.3 


1.4 


1.4 


1.4 


1.4 


1.4 


3.0 


1.70 


4.0 


4.5 


4.8 


4.9 


4.9 


4.9 


4.9 


2.8 


1.70 


3.6 


4.4 


5.4 


5.8 


5.9 


5.9 


5.9 


2.8 


1.75 


5.2 


6.6 


7.1 


7.7 


7.8 


7.8 


7.8 


2.7 


1.65 


3.5 


4.5 


5.0 


6.1 


6.2 


6.2 


6.2 


2.6 


1.65 


3.2 


4.5 


5.2 


6.0 


6.1 


6.1 


6.0 


2.5 


1.65 


3.2 


4.4 


6.8 


7.5 


7.6 


7.6 


7.6 


2.4 


1.72 


4.5 


7.2 


12. 


16. 


17. 


17. 


17. 


2.2 


1.70 


6.2 


15. 


20. 


43. 


53. 


54. 


80. 


2.0 


1.60 


5.0 


13. 


54. 


139. 


199. 


211. 


213. 


0.0 


1.45 


28. 


36. 


119. 


524. 


1156. 


1505. 


1570. 



Table 4: Dependence of relaxation times r on Am^ in the ideal MG algorithm with lexicographic 
SOR for computing propagators of staggered fermions on 12^ lattices. (The values given for 
P — 2.5 were obtained in the same gauge field configuration.) 











T 


for Am 


2 _ 








cu 


1 


10-1 


10-2 


10-3 


10-^ 


10-5 


10-6 


oo 


1.32 


1.5 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


3.0 


1.65 


2.9 


3.2 


3.3 


3.5 


3.5 


3.5 


3.5 


2.7 


1.65 


2.9 


3.4 


4.0 


4.5 


4.5 


4.5 


4.5 


2.5 


1.72 


4.2 


4.9 


6.9 


8.4 


8.7 


8.8 


8.8 


2.5 


1.65 


3.2 


3.8 


6.9 


8.8 


9.3 


9.4 


9.4 


2.4 


1.65 


3.3 


5.0 


15.4 


23.6 


27.7 


28.0 


28.0 
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